(3+l)-TQFTS AND TOPOLOGICAL INSULATORS 

KEVIN WALKER AND ZHENGHAN WANG 

Abstract. Levin- Wen models are microscopic spin models for topological phases 
of matter in (2 + l)-dimension. We introduce a generalization of such models 
to (3 + l)-dimension based on unitary braided fusion categories, also known 
as unitary premodular categories. We discuss the ground state degeneracy on 
3-manifolds and statistics of excitations which include both points and defect 
loops. Potential connections with recently proposed fractional topological insu- 
lators and projective ribbon permutation statistics are described. 

Key words: TQFT, topological insulator, premodular category 

Pacs: 02.40.Pc, 03.65.Vf, 03.75.Lm 



1. Introduction 

In a remarkable paper, M. Levin and X.-G. Wen defined a family of rigorously 
solvable lattice spin Hamiltonians for a large class of 2D topological phases of 
matter based on string-net, called the Levin- Wen model [LW] . The Levin- Wen 
model takes a unitary fusion category C as input, and as output, realizes the 
topological phase effectively described by the topological quantum field theory 
(TQFT) based on the Drinfeld center Z(C) of C. In Section V of their paper, Levin- 
Wen model is generalized to 3D (and > 4 dimensions) using unitary symmetric 
fusion categories. The 3-D Levin- Wen model realizes all discrete gauge theories 
coupled with bosons and fermions. In this paper, we show that a generalized 
string-net model exists for all unitary braided fusion categories, which include 
the unitary symmetric fusion categories and unitary modular categories as two 
special cases. The case for unitary modular categories is a generalization of the 
(3 + 1)-BF theories. A unitary braided fusion category, also called a unitary 
premodular category, is some non-trivial product of a discrete gauge theory with a 
unitary modular category. Since unitary modular categories are algebraic theories 
of anyons, therefore our new models can be thought as discrete gauge theories 
coupled with anyons. 

As a generalization of the Levin- Wen model, our Hamiltonians also consist of two 
kinds of commuting projectors, and are stable under small, yet arbitrary pertur- 
bations. On the 3-sphere, the ground state manifold is non-degenerate. In general 
the discrete gauge theory part of a unitary premodular category corresponds to a 



The second author is grateful to Xiao-Gang Wen for very insightful comments on an earlier 
draft of the paper. We thank L. Chang for drawing most of the pictures. 

1 



2 



KEVIN WALKER AND ZHENGHAN WANG 



finite abelian group G, and we expect the ground state manifold of our Hamilton- 
ian on a general 3-manifold X is isomorphic to C[Hi(X;G)]. The ground state 
manifold on any 3-manifold X are examples of 3-D-error correction codes. It will 
be interesting to know if there are examples of self-correcting quantum memories 
for some theories on certain 3-manifolds. Pointed excitations in our models are 
still bosons and fermions. What interesting is the existence of extended excitations 
such as loops and ^-graphs. The mutual statistics of pointed excitations and the 
extended excitations are more general than bosonic and fermionic statistics. 

When the time reversal symmetry of fractional topological insulators is bro- 
ken, fractional topological insulators can be connected to some topological phases 
including the trivial one. Therefore, 3D fractional topological insulators can be 
considered as topological phases with symmetry [Wenlj . We conjecture that our 
(3 + l)-TQFTs are the underlying topological orders for 3D topological insulators 
and their generalizations. A classification of the compatible symmetries in our 
models using the projective symmetry group [Wen2] should reveal a connection 
to topological orders with symmetry. A classification of time-reversal symmetries 
within a given unitary premodular category would lead to a classification of all 
3-D-topological insulators with the same underlying topological order. 

2. Topology in Condensed Matter Physics 

As a branch of pure mathematics, topology is the study of spaces regardless of 
metric. Therefore, topology in physics seems to be a strange occurrence because 
distance and time are of paramount importance in physical measurements. But 
the first application of topology in physics predates topology and goes back at 
least to 1833 when Gauss revealed his beautiful formula for the linking number. It 
was argued that the formula originated first from Gauss's study of the tracking of 
asteroids and comets, and later he applied his formula to electromagnetism |Epp] . 

Suppose L\, L 2 are two disjoint simple closed curves in M 3 . If there is an electric 
current of strength j 1 in the wire Li = {x'(s)\s G S 1 }, then it generates a magnetic 
field 

= W 1 I (x'-x)xdx'^ 
4tt J Ll \x'-x\ 

by the Biot-Savant law for x G IR 3 \Li. By the Maxwell equation § L B(x)dx = fioj 
for L 2 = {x(t)\t G S 1 }, it follows that 

/ioji / f (x' — x) • {dx' x dx) 



Air 



l 2 J Li K ~ ^\ 

Therefore, an electromagnetic definition of the linking number between L\ and 
L 2 will be the ratio 4^. Gauss's formula has a geometric explanation as follows. 

Taking two points x',x on L\,L 2 and normalizing the line segment x'-x to ifrEfj ? 

then we obtain a map from the abstract torus S 1 x S 1 of s, t parameters to the unit 
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sphere 5 2 C R 3 . Gauss's formula is simply the degree of this map. Though linking 
number is a very useful invariant for linkage, Maxwell, who were very interested in 
knot theory, seemed to be the first to find two un-separable disjoint simple closed 
curves with linking number equal to |Epp| . 

An important topic in topology is the study of special spaces called manifolds. 
An n-manifold M is a space that is locally the Euclidean space M. n up to homeo- 
morphism. Simple examples are n-spheres: 1-sphere is the circle, and 2-sphere is 
our ordinary sphere. It is a challenge to visualize manifolds beyond dimension 2, 
and see their relevance to condensed matter physics. After all, all physical experi- 
ments are carried out in our 3-space 1R 3 , and most 3- manifolds cannot be embeded 
into M 3 . But complicated manifolds could arise in condensed matter physics at 
least in two different situations: as subsets of R 3 with complicated boundary con- 
ditions, or as configuration spaces. Just as every orientable surface is a polygon in 
the plane with glued sides, every orientable 3-manifold can be obtained by iden- 
tifying pairs of faces of a polyhedron (solid). A familiar example is the 3-torus 
as a cube with periodic boundary identifications. Another famous example is 
the Poincare homology 3-sphere obtained by identifying faces of a dodecahedron. 
Topologists use topological invariants to distinguish manifolds. Topological invari- 
ants in physics often arise classically by integrating local geometric quantities such 
as curvature, or quantum mechanically by path integrals of total derivative terms, 
usually dropped, in the action functional. 

The discovery of the fractional quantum Hall (FQH) liquids, and recently the 
topological insulators stirred great interests in topology in condensed matter physics. 
Witten-Chern-Simons theories as effective theories for FQH liquids and Chern 
numbers used in the study of topological insulators represent two kinds of topo- 
logical invariants: quantum and classical. While there is no clear-cut separation, 
by quantum invariants we mean invariants of spaces obtained as path integrals of 
some TQFTs. Famous examples are the Jones polynomials at roots of unity from 
Witten-Chern-Simons theories. Classical topological invariants include homotopy 
groups, (generalized) (co) homology groups such as K-theory, and characteristic 
classes such as Chern classes. Using this rough division, the Jones polynomial at 
4-th root of unity appeared for the 5 / 2-fractional quantum Hall effect is a quantum 
invariant, while the Chern number appeared in topological insulators is a classi- 
cal invariant. But some invariants defined quantum mechanically turn out to be 
classically determined. 

One of the frontiers in topological phases of matter is the understanding of in- 
teracting 3.D-topological insulators [MQKS| [SBMSJ. Interactions can be thought 
as dynamical entanglement. When entanglement becomes long-rang, topological 
symmetry could emerge, and an effective description by a TQFT is possible. There- 
fore materialized topological symmetry can be described by the tensor category 
encoding the effective TQFT. Non-interacting topological insulators are believed 
to be modeled by (3 + 1)-BF TQFTs [CM], whose path integrals are determined 
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by classical topological invariants. We ask: how many (3 + l)-TQFTs do we know 
and which might be related to interacting 3-D-topological insulators? 



3. (3+l)-TQFTs 

Dimension 4 is different. The Euclidean space M. is the only Euclidean space 
that has more than one smooth structure. Roughly, this means that there is more 
than one way to do calculus on IR 4 . And not just one more — there are infinitely 
many more different ways to do calculus on M 4 [Gomj . The classification of smooth 
4-manifolds is one of the most difficult problems in mathematics (there is another 
flavor of 4-dimensional topology: the classification of topological 4-manifolds. Due 
to M. Freedman's work, the landscape of topological 4-manifolds is more or less 
understood FQ|.) Therefore, it is not surprising that there are not many smooth 



topological invariants, classical or quantum, in dimension 4 that detect smooth 
structures. In particular, (3 + I)-TQFTs are very difficult to find because they 
could generate smooth topological invariants of 4-manifolds. The most powerful 
one for AD topology is the first TQFT invented by E. Witten in 1988: an N = 2 
supersymmetric Yang-Mills theory which reproduces the Donaldson invariant of 
4-manifolds |Wittlj . 

As a first approximation, TQFTs are quantum field theories such that the path 
integrals Z(M) = Tr(e ltH ) for a theory with Hamiltonian if on a space manifold N 
is a topological invariant, where M = N x S 1 is the space-time with periodic time. 
In particular Z(M) should be independent of time t. (2 + l)-Witten-Chern-Simons 
type TQFTs are those with H = 0, but the (3 + l)-Witten-Donaldson TQFT is 
more subtle. The trace is interpreted as a super-trace Tr((— l) F e' ltH ), where (— 1) F 
is the fermionic parity operator, so time dependence could be canceled. Moreover, 
Donaldson invariant is defined only for space-time manifolds satisfying certain 
topological restrictions [DKJ. 

Two interesting families (3 + l)-TQFTs are discrete gauge theories, and BF 
theories [DWj[BaezJ. Both families of TQFTs give rise to topological invariants 
that are determined by classical homotopy invariants. The same should be true 
for all the (3 + l)-TQFTs based on unitary braided fusion categories. 

Witten-Donaldson theory is not determined by homotopy invariants, and detects 
smooth structures of 4-manifolds. The Witten-Donaldson TQFT is actually a 
partial TQFT because it is defined only for 4-manifolds with b£ > 1. For example, 
it is not defined for the 4-sphere S 4 . Very recently, Witten defined another (3 + 1)- 
TQFT: an iV = 4 supersymmetric Yang-Mills theory |Wit t2j. It is not known if 
this TQFT can detect smooth structures of 4-manifolds. There are some other 
proposed (3 + l)-TQFTs or partial TQFTs, but not much is known about them 
for the detection of smooth structures. 
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Physicists describe fermions using Grassmann numbers, which are differential 
1-forms, in the second quantized framework. By De Rham theory, closed differ- 
ential forms represent cohomology classes, thus fermions are aware of the topol- 
ogy of the space that they occupy. Therefore, topological quantum field theo- 
ries with fermions are potentially very different from purely bosonic ones such 
as (2 + l)-Witten-Chern-Simons theories (Reshetikhin-Turaev TQFTs mathemat- 
ically). Witten-Chern-Simons theories are bosonic in the sense that the path 
integrals contain no Grassmann numbers. 

3.1. Discrete Gauge Theories. Discrete gauge theories based on finite groups 
G are the best understood examples of higher dimensional TQFTs [Freed] . For a 
(3 + 1)-TQFT based on G, the path integral of a 4-manifold W counts the number 
of representations from the fundamental group of W to G. The Hilbert space as- 
sociated to a 3-manifold X is spanned by the conjugacy classes of homomorphisms 
from the fundamental group of X to G. 

3.2. BF Theories. Given a Lie group G, and a principle G-bundle P over a 
4-manifold W, the (3 + 1)-BF theory is a TQFT based on the action 



where B is a adP-valued 2-form, F the curvature of a connection A on P, and 
A is a coupling constant |Baezj . When G = GL(4,C), A plays the role of the 
cosmological constant. The path integral for a 4-manifold W is e l ^ <T<yW \ where f3 
is some constant independent of W, and a(W) is the signature of W. The Hilbert 
space associated to a 3-manifold X is always 1-dimensional. 

3.3. Cohomological Field Theories. Witten-Donaldson type (n + 1)-TQFT 
are cohomological in the sense that the vector space associated to an n-manifold 
M is related to the cohomology of some moduli space. For Donaldson theory, 
fix a compact Lie group G and a principle G-bundle P over a 4-manifold W. 
Then the moduli space is the space of anti-self dual instantons, i.e., solution to 
the Yang-Mills equation F A = for connections A on P. Donaldson invariant 
roughly counts algebraically the number of solutions. An important ingredient in 
the formulation of Donaldson theory is supersymmetry (topological twist, which 
is also very important, will not be discussed here.) Witten formulated Donaldson 
invariant as correlation functions of an A^ = 2 supersymmetric Yang-Mills theory 
[Wittlj . A beautiful illustration of the role of supersymmetry is the following 
formulaQ to count algebraically the number A/ of solutions to f(x) = of a 
generic smooth function / : [0, 2] — >■ [—1, 1] such that /(0) = — 1, /(2) = 1: 



The second author learned the formula from J. Maciejko who attributed it to Witten. 
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where \ and if) are Grassmann variables. 

An important event in 4D topology is the discovery of Seiberg-Witten equation 
around 1994 |Witt3j . Seiberg-Witten invariant is conjectured to give the same 
topological information about 4- manifolds as Donaldson invariant, and is much 
easier to work with mathematically Seiberg-Witten theory is the infra-red limit 
of the N = 2 supersymmetric Yang-Mills theory. Given a spin c -structure on the 
4-manifold W, there associated two spinor bundles A ± and a line bundle L over 
W. The Seiberg-Witten equation is for a pair [A, if}), where A is connection on L 
and ip a positive spinor: 

d a V = o,f! = -<7(v,^), 

where Da is the Dirac operator, and a is some paring for positive spinors. The 
Seiberg-Witten invariant basically counts the number of solutions. 

4. Discrete Gauge theory coupled with anyons 

4.1. Rigorously Solvable Hamiltoninans. The existence of rigorously solvable 
Hamiltonians follows from a general mathematical theory of picture TQFTs |Wal] . 
Such (3 + l)-TQFTs based on unitary braided fusion categories will be studied 
mathematically in [WW] . In this paper, we will write down explicit Hamiltonians 
for those (3 + l)-TQFTs and quote necessary results from [Wal] [WW] . 

For convenience, we will describe the model on the cubic lattice in detail, and 
only sketch the model for general 3-manifolds. 

4.1.1. Algebraic data. There are several equivalent definitions of a unitary braided 
fusion category (UBFC). We refer the interested readers to [Wanj . While ev- 
ery definition is complicated, there is a graphical calculus and interpretation for 
a UBFC as an algebraic theory of anyons which make the axioms reasonable. 
The graphs involved can be thought as anyon trajectories and interesting topo- 
logical changes correspond to physical events such as creation/annihilation, fu- 
sion/spliting, braiding and twisting. Given an anyoic system with particle types 
a, b, c, ... which form a label set L. The number of particle types is finite and 
called the rank of the UBFC. A particular convenient way to present a UBFC 
is through three sets of numbers: {t a = ±l} a( =L, {F^ m }, a, b, c, d, m, n G L, and 
{R® b },a,b, c G L. The numbers {t a = ±l} ag £ encode the creation/annihilation 
structures. The so-called 6j symbols {F^ m },a,b,c,d,,m,n G L are solutions to 
the pentagon equations, while {R^}, a,b,c G L encode the braidings which satisfy 
the hexagon equations. The three set of numbers are not independent as the cre- 
ation/annihilation, fusion/spliting, and braiding have to be compatible. (Strictly 
speaking the 6j-symbols are lOj symbols and the braidings R^ b are unitary matri- 
ces when some fusion multiplicity is > 1. For simplicity, we assume all UBFCs are 
multiplicity- free.) As a consequence, any planar trivalent graph with braidings can 
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be evaluated to a number, which is the topological amplitude of the represented 
physical process. 

The 6j-symbols {F°[ h r c im } can be organized into matrices, called F-matrices, by 
the following diagram: 

a b c a b c 



= E 



rpabc 
n d;nm 



n 



d d 

Similarly, the braiding eigenvalues are defined by the following diagram: 

a b 



R. 



ab 



C C 

Besides multiplicity-free in the fusion rules, we will also assume that all labels 
are self-dual, so edges in our lattices are not oriented. Furthermore, we assume 
edges around any vertex can be bent as long as we do not introduce any crossings. 
Examples of such theories are premodular categories from the Temperley-Lieb alge- 
bras [WanJ. The Hamiltonian below generalizes to the general case with adequate 
notation. 



4.1.2. Cubic lathee model. Given {F^ m }, {Rf} of a UBFC with label set L 
(strictly speaking, we should choose a set of representative simple objects). Let 
C L be the Hilbert space spanned by all labels. Just as in the Levin- Wen model, 
it is convenient to work with trivalent graphs, therefore we first resolve the cubic 
lattice C into a trivalent lattice. There are many ways to do it, and the resulting 
theories are all equivalent as each one is equivalent to the continuous limit. At 
each 6-valent vertex, we resolve it as follows: 





As a result, each 6-valent vertex of the cubic lattice C is split into four trivalent 
ones with a Z 3 -symmetry. We assume periodic boundary conditions, so our lattice 
is in the 3-torus. 
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The Hilbert space V of the model is spanned by all labelings of the edges in 
the resolved cubic lattice, denoted by T C - Equivalently we assign a qudit C L to 
each edge of Tc including the new ones. There are two kinds of terms in the 
Hamiltonian: the vertex type and the plaquette type. In order to write down the 
Hamiltonian explicitly, we need to project our lattice onto a plane. We look at 
the lattice along the (1,1, Indirection, and therefore the resolved cubic lattice T c 
is projected to a plane perpendicular to the (1, 1, Indirection. The microscopic 
Hilbert space is V = <S>eerc^ L - The Hamiltonian H is of the form —J\ ^2 veTc A v — 
J2 X] P ec -^P' w here v ranges over all vertices of Tc including the new ones, p 
ranges over C sq — all plaquettes which correspond to the original squares of the 
cubic lattice C, and J±, J 2 are coupling constants. 

The microscopic Hilbert space V = ® e erC L has a natural basis: the 1-skeleton 
of Fq with a label on each edge. Then for a vertex v of Tc, the term A v acts 
a basis as follows: A„|\I/) = if the three labels around v obey the fusion 
rules; otherwise A v \^/) = 0. A plaquette term is much more complicated, we 
explain below how it is derived and give its formula in terms of 6j-symbols and 
braiding eigenvalues. 

Our resolution of the cubic lattice to a trivalent one is very symmetric. There 
are three kinds of plaquettes in C sq : those corresponding to squares in the planes 
parallel to x-y, y-z and x-z plane. Below are the pictures of three representatives 
of such plaquettes. 




y 



Note that each plaquette p is a decagon (although we will call such a plaquette a 
decagon, it is not always planar. It is a surface bounded by a polygonal path with 
10 edges: 4 of them correspond to the 4 sides of the original square of the cubic 
lattice and 6 new ones from the resolution.) There are also 10 adjacent edges. 
The plaquette term indexed by p is obtained by placing a simple loop labeled by 
a projector uq onto the plaquette along the edges as below. 
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Physically the projector ooo enforces the total flux through p to be a transparent 
label (see Section H~5l below for the definition). Such a projector is formally written 
as Co> = J2s€L W s i w h ere d s is the quantum dimension of the label s and D 2 = 
J2seL Adding such a loop with the projector u will not change the topological 
amplitude of a basis |^}. This can be seen by expanding ui into J2sgl W s anc ^ 
noticing that a contractible loop labeled by s is evaluated to d s . A formula for a 
plaquette term is then obtained by evaluating the same projector in a different way 
using 6j symbols and braiding eigenvalues. For analogous derivations of similar 
terms in Levin- Wen model, see page 100 of [WanJ. 

Due to the regularity of the cubic lattice and the symmetry of our resolution, 
we need only to write down one plaquette term. We choose to write the formula 
for the plaquette in the x-y-plane, denoted as p xy . We could equally work with 
the one in the x-z-plane p xz or the one in the y-z-p\aiae p yz . 

To write down such a formula, we denote the basis element that labels the 10 
edges of p xy by abcdpqruvw and their 10 adjacent edges by a'b'c'd'p'q'r'u'v'w' as 
in the following picture by abcdpqruvw)- Labels of edges that are not named 
remain the same in all computations. Our convention is that the edge not on the 
decagon p , but next to the edge of the decagon labeled by I, is labeled by V . 



The plaquette term B p will map the basis vector \^ Px , a bcdpqruvw) into a big 
linear combination of basis elements, where the labels a,b,c,d,p,q,r,u,v,w are 
replaced by new labels a", b", c", d", p", q", r", u", v", w". In the following, 




a in a in a a a a a 

a ,o ,c , a ,p ,q ,r ,u ,v ,w 



n 
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will be abbreviated as a", w" . Therefore, all we need are the coefficients „, 



ill Bp | p xy ,abcdpqruvw) 



" in" 



d. 



1)2 Bp xyt a",...,w"\^p xy ,a",...,w") 
seL a",...,w"eL 



Recall that uo is the formal sum J2 seL jfzs. The operator B p = J2 seL j^B p 
is a sum of operators , where -B* is the operator that corresponds to the 

rxy rxy 

simple loop labeled by s. Hence it suffices to know the coefficients B 1 ! in 

B s p jy Pxy , abcdpq ruv W ) = Ea'',...,«» B J J „a'',...y'l^ I a'' 1 ...yl In terms of 6 J symbols 
and braiding eigenvalues, we claim 



n p xy ,a",...,w" 



rjq'b TDc'r Dl'°" r>c'r" rpa"sp rpp"sq -pq" sb rpb" sc rpc" sr rpr" su rpW'sd rpd" sv rpv" sw rpW'sa 
n q n c n q „ tt d , r a /. ap /i-T p '- pq " p q '- q b" r b';bc" r c';cr" r r';ru" r u';ud" P d';dv" r v';vw" r w';wa"- 

To derive this formula, we first twist the labeled graph representing the basis 
\^p x ,abcdxyzuvw) around the two vertical edges as below. 



su jpu"sd Tpd" sv rpv" sw rpw"sa 



v' w 




This multiplies \^ Pt abcdxyzuvw) by R q 6 -R"' c . Then we fuse the simple loop labeled 
by s with the edge labeled by a as shown below. 




v' w \ a a" a/. V p' 

Next a sequence of F-moves brings the s-labeled strand counter-clock-wise along 
the boundary of the decagon p through all the trivalent vertices one by one. Each 
time when the s-labeled strand passes a trivalent vertex on p, an F-move is used. 
Due to the two introduced twists, we do not need to use braidings when we perform 
all the F-moves. 
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Finally, the s-labeled strand returns to the edge labeled by a with a bubble. 




Removing the bubble and twisting back the two vertical edges, we arrive at the 
formula above. 

4.1.3. Dual cellulations . The cubic lattice is so regular that a lot of technical 
difficulties disappeared. Since there is substantial topology involved for the general 
case, we are content with a sketch of the procedures. For the basic topology 
involved, see the book |BP] . 

Trivalent graphs are generic in surfaces under perturbations. Genetic 2-dimensional 
polyhedra are simple polyhedra which are generalizations of trivalent graphs. For 
a general 3-manifold X, our Hamiltonian is given on a branched standard spine 
of X. An elementary way to present a 3-manifold X is by a triangulation of X. 
A triangulation A of a 3-manifold X is a collection of tetrahedra {Aj} such that 
tetrahedra are glued together along their faces. Our model is conveniently defined 
using the dual cellulation Ta of the triangulation A. In two dimension, recall 
that in the dual cellulation of a triangulation of a surface, each triangle becomes a 
vertex, an edge still an edge, while a vertex becomes a cell (polygon). Notice that 
the 1-skeleton of the dual cellulation of a surface is always a trivalent graph. In 
3D, a vertex of the dual cellulation Ta is the center of a tetrahedron. Two vertices 
of Ta are connected by an edge if their corresponding tetrahedra share a face. A 
face of Ta is dual to an edge of A, and a 3-cell (solid) dual to a vertex of A. The 
dual cellutation of a 3-manifold is an example of a simple polyhedron. Our model 
can be easily generalized to any branched standard spine of a 3-manifold X, which 
is much more convenient to work with in practice. 
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In order to define our Hamiltonian, we need an extra structure, called an oriented 
branching, on the triangulation A, or equivalently on i~A- We will define it for 
the triangulation A, but the translation to Ia is straightforward. An oriented 
branching on a triangulation A of X is an assignment of an arrow to each edge of 
A so that in the three edges of each triangle of A there are exactly two consistent 
arrows, i.e., the three arrows of any triangle never form a cycle. Every oriented 
3-manifold A has a triangulation with an oriented branching |BP] . An oriented 
branching uniquely determines an ordering of the 4 vertices of each tetrahedron if 
the arrows go from lower numbered vertices to higher numbered ones. 

Using an oriented branching, we can identify each tetrahedron in a triangulation 
A of A with the standard tetrahedron in M 3 . The 1-skeleton of the dual cellulation 
is a 4-valent graph. We can resolve this 4-valent graph into a trivalent graph in 
a standard way. The Hilbert space of our model is again the tensor product of 
all qudits C L over all edges of the resolved 1-skeleton of the dual cellulation IA 
including the new ones. There are again two kinds of terms: vertex and plaquette 
types. The Hamiltonian is H = — J2 v er ^ v ~ J2 p er B p , where v ranges over all 
vertices including the new ones, and p ranges over all 2-cells of the dual cellulation 
(plaquettes). The oriented branching allows us to have standard local models of 
the standard spine in M 3 , therefore we can use a similar procedure as for the cubic 
lattice to write down the plaquette terms. 

4.2. Ground State Manifold. Given an oriented 3-manifold A, consider the 
infinite dimensional vector space A(X) generated by all colored string- nets in A. 
Let A(X) be the quotient space of A(X) by all local relations derived from the 
UBFC. (Strictly speaking, the string-net strands are ribbons and the vertices are 
rigid. We refer the interested readers to [WW J for a mathematical discussion.) 
Then A(X), called the skein spaces, are isomorphic to the ground state manifolds 
of the spin Hamiltonians, which are examples of error-correction codes. 

4.3. Statistics of Excitations. Elementary excitations in the spin models in- 
clude pointed particles, loop-defects, ^-defects, and more general defects. Their 
types and statistics can be described using representations of certain cylindrical 
categories related to boundary conditions for 3-manifolds. 

We consider two kinds of boundary conditions for a 3-manifold A with a bound- 
ary surface Y: crude and topological. A crude boundary condition c on Y is a 
finite collection of points labeled by objects of C. Let A(X, Y; c) be the vector 
space generated by all colored string-nets in A which terminate at c. Then the 
relative skein space A(X, Y; c) is the quotient of A(X, Y; c) by local relations. We 
define a category A(Y) for each surface Y (including the empty one) as follows. 
An object of A(Y) is a crude boundary condition c. The morphism space from a 
boundary condition c\ to C2 is the relative skein space A(Y x [0, 1]; ci, C2). There- 
fore, a morphism from an object C\ to another object c<i is represented by a linear 
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combination of string-nets in Y x [0, 1] that terminate at c\ and C2 in Y x and 
Fx 1, respectively. 

Just as algebras, the linear categories A(Y) have representations. Any repre- 
sentation of A{Y) is called a topological boundary condition. An irreducible rep- 
resentation corresponds to an elementary excitation whose boundary is Y . The 
statistics of elementary excitations can be computed by using the relative skein 
space with topological boundary conditions. (2 + l)-dimensional analogues are 
discussed in [FNWWj . 

4.4. Modular Category. The SU (2)-Witten-Chern-Simons theories can be pro- 
moted to (3 + l)-TQFTs, called Crane- Yetter TQFTs [CV]. The Crane- Yetter 
(3 + l)-TQFTs and their generalizations based on premodular categories induce 
representations of the motion groups — generalizations of the braid group to more 
general mapping class groups. 

The generalization of the braid group to extended objects in 3D is conceptually 
straightforward. The generalization to small loops is called the loop braid group 
[Lin] . Consider finitely many small 3-balls {Bf} inside the standard 3-ball B 3 . 
Fix the equator Si of each small 3-ball Bf. Then the mapping class group of 
5 3 \{S'j} — the self-diffeomorphisms of -B 3 \{S'j} fixing the outside boundary of B 3 
modulo deformations, is the loop braid group. Given a group, it is not always easy 
to find a presentation, i.e., a set of generators and relations. For the loop braid 
group, this is done in [Lin]. 

Suppose (V, Z) is a (3 + 1)-TQFT based on a modular category C. For an 
oriented closed 4-manifold W, the topological invariant Z(W) is e l ^ ca ^ w \ where 
c is the topological central charge of C, and a{W) is the signature of W . For 
an oriented closed 3-manifold, the vector space is always 1-dimensional. When an 
oriented 3-manifold X has a boundary Y, then the topological boudnary condition 
for Y is always trivial. The vector space associated to X with this trivial boundary 
condition on Y is isomorphic to Vrt(Y) — the vector space associated to Y in the 
Reshetikhin-Turaev (2 + 1)-TQFT based on C. 

Mathematically, the category A(Y) is trivial up to Morita equivalence |WWj . 
Consequently for any oriented closed surface Y, there are neither non-trivial par- 
ticle excitations nor non-trivial excitations of extended objects. This is analogous 
to the situation of (2 + l)-TQFTs based on modular categories such as (E s )i- For 
modular categories from Chern-Simons theories, our models gap out the F A F- 
theories in the bulk, while the Chern-Simosn theories on the boundary survive. 

4.5. Mathematical Underpinning and Continuous Limit. The conceptual 
underpinning of the 2D Levin- Wen model is two mathematical theorems: The 
quantum double Z{C) of a unitary fusion category C is always modular, and the 
Turaev-Viro TQFT based on C is equivalent to the Reshetikhin-Turaev TQFT 
based on Z(C) [Muej jTVi] |BenK] . Therefore, Levin- Wen model is a Hamiltonian 
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realization of both theorems simultaneously. The mathematical theory behind the 
new models should be a generalization of these theorems to (3 + l)-dimension. 

There is a structure theorem for premodular categories. An object x in a pre- 
modular category is transparent by definition if s xy = d x d y for any object y, where 
s xy is the topological invariant of the Hopf link colored by x and y. Transparent 
elements of a premodular category C form a symmetric fusion category Sc- By a 
theorem of Deligne, every symmetric fusion category is equivalent to the represen- 
tation category of a pair (G, /i), where G is a finite group and /i is a central element 
of G of order < 2 (see [Ostj ). When /i is the identity, then there is a quotient Qc 
of C by Sq, which is modular. When /i is not the identity, a generalized quotient 
exists. Therefore, essentially C is some kind of extension of the quotient Qc by Sc- 

The spin lattice models have continuous limits which are TQFTs based on pic- 
tures. A framework for constructing picture TQFTs is formulated in |Walj . The 
ribbon graphs based on a premodular category are examples of a system of fields 
in the sense of |Walj and hence lead to a (3 + l)-TQFT. In general, picture TQFTs 
have state-sum formulation, which can be realized by spin models. 



4.6. Holographic Resolution of Anomaly. (2+l)-Witten-Chern-Simons (WCS) 
theories have anomaly in the sense that the path integral for closed 3-manifolds 
are not well-defined unless the 3-manifolds are endowed with some extra structures 
such as 2- framing [Witt4] . In dimension 2 the anomaly is manifested in the chiral 
central charge of the boundary CFT. The (3 + 1)-TQFT based on a modular cate- 
gory can be regarded as a holographic resolution of the anomaly: the (2 + 1)-WCS 
theories are really (3 + l)-TQFTs. Our model leads to a holomorphic tensor net- 
work representation for Chern-Simons theories, which is probably related to some 
generalization of the AdS / CFT correspondence. 

Given a modular category C, the Reshetikhin-Turaev (RT) TQFT leads to a 
topological invariant of a 3-manifold X with some extra-structure. A convenient 
way to encode this extra structure is by an integer n. Therefore, by an extended 
3-manifold, we mean a pair (X, n), where X is an oriented 3-manifold and n an 
integer. For a surface Y, the extra structure can be given a Lagrangian subspace 
L in Hi(Y;Q). RT TQFT will associate a vector space to the pair (Y, L). From 
the (3 + l)-point of view, the extra structures n and L give instructions to finding 
manifolds one-dimensional higher. Then the associated path integral and vector 
space for these one-dimension higher manifolds are exactly the path integral and 
vector space from RT TQFTs. Specifically, given an extended 3-manifold (X, n), 
choose an oriented 4- manifold W bounding X and a(W) = n, then Z^ + i(W) = 
Zrt{{X,ti)). For an extended surface (Y,L), choose a 3-manifold X so that the 
kernel of the inclusion of H X {Y] Q) in H X (X] Q) is L, then Z 3+1 (X) = Z RT ((Y, L)). 
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5. Applications 

5.1. 3D Topological Insulators and BF theories. There are several proposals 
for the effective TQFTs that model topological insulators |QHS|[CM] . The mod- 
eling with BF theories is of particular interests to us because (3 + 1)-BF theories 
are related to (3 + l)-TQFTs based on unitary modular categories. 

Consider the (3 + 1)-BF theory with G = SU(2). By integrating out the _B-field 
when A 7^ 0, we have 

Z(W)= [ [ e l fwTr(BAF+±FAF) DBDA= ^ f e -%f w Tr(FAF) DA 

J Jb,a vA J a 

Comparing with the S77(2)-WCS theory, we see that ^ corresponds to the level 
k in S77(2)-WCS theories. It is generally believed that the path integrals in 
the top dimension of a TQFT determine the extended TQFTs down to at least 
codimension= 2. If so, then (3 + 1)-BF theories are the same as WCS theories 
promoted to (3 + 1). When A — > 0, which is the same as k — > oo, the limit theory 
is a semi-classical one. 

As is explained in earlier sections, it is possible to compute statistics of the 
elementary extended excitations in the (3 + l)-TQFTs, though in practice it is not 
an easy task. Physically, we can couple the rank= 2 tensor field in equ. (30) of 
[CM] with extended objects such as a line field for G = SU (2) (choose a A corre- 
sponding to some level k, say k = 2, in S77(2)-WCS theory.) Conjecturally, the 
resulting statistics will be the same as computed from the mathematical theories. 
It will be interesting to check some examples. 

More interestingly, there are several proposals for potential fractional topologi- 
cal insulators |MQKS| |SBMS] . As alluded in the introduction, we believe that the 
underlying topological orders in the fractional topological insulators when symme- 
try is broken will correspond to some (3 + l)-TQFTs based on UBFCs. It will 
be interesting to find out what are the corresponding UBFCs for the proposed 
fractional topological insulators. 

5.2. Projective Ribbon Permutation Statistics. In |F HNQWW| , the possi- 
bility of ribbon permutation statistics is studied for a collection of confined pointed 
excitations, called hedgehogs. If the hedgehogs become deconfined in some related 
theory, then an effective description by a (3 + 1)-TQFT is a possibility. We con- 
jecture that one such possible deconfined phase could be a (3 + 1)-TQFT based 
on the following unitary braided fusion category. 

Fermionic quantum Hall liquids can be described using spin modular categories 
[RWJ. For |-FQH liquids, the fermionic Moore- Read state with 6- fold degeneracy 
on the torus is covered by a rank= 12 UBFC. This category is the even sector 
of I sing x Z§: the direct product of the Ising theory and a modular theory with 
fusion rules Zg. The Ising theory has anyon types {l,a,ip} and the anyon types 
of the modular Z 8 theory are denoted by {0, 1, 7}. The anyon / = ip <g> 4 is a 
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fermion. The even sector, the anyons which are local with respect to the fermion 
/, consists of {1 (g) i, ip <g) i} for i =even and {a (g> i} for i =odd. This even sector 
is a rank= 12 UBFC, therefore it leads to a (3 + 1)-TQFT. 



6. Further discussion 

The most general TQFTs are given by universal manifold pairing [FKNS W Wj . 
This is tautology: the invariant of a manifold is itself considered as a vector in 
some vector space. Dimension 4 is again different from lower dimensions because 
the pairing has null (or light-like) vectors, while there are no such null vectors in 
dimensions 1, 2, 3. If we restrict to unitary and anomaly-free TQFTs, then some 
smooth structures of some 4-manifolds cannot be distinguished by the universal 
manifold pairing in dimension 4. 

We expect a true 3D generalization of the Levin- Wen model will be based on 
a generalization of unitary fusion category to a unitary fusion 3-category and 
membranes. Thus it might seem strange that interesting (3 + l)-TQFTs can be 
constructed based on string-nets. But the string-net that we are using are mem- 
branes: the strings are ribbons, so we are working with very simple membranes. 
Still why the 6j symbols are useful algebraic input because the algebraic input 
should be solutions to a generalization of the pentagons to 3D? Braided tensor 
categories are examples of 3-categories. A curious fact is that the pentagons wear 
two hats: in 2D, it is the algebraic equation for diagonal flip, while in 3-D, it is 
the Pachner 2-3 move. Presumably this is a manifestation of the relation between 
2D CFTs and 3D TQFTs. 

If no fermions are involved, it is not clear if a high category generalization of 
Levin- Wen model would produce (3 + l)-TQFTs to distinguish smooth structures 
of 4-manifolds. Based on universal manifold pairing and Witten-Donadlson theory, 
we would speculate that fermions are important for formulating (3 + l)-TQFTs 
capable of distinguishing smooth structures. In 2D, a generalization of Levin- Wen 
model to include fermions is proposed in [G WWj . which shows that topological 
orders in systems with fermions are strictly richer than purely bosonic systems. 

The N = 2 supersymemtric Yang-Mills theory for Witten-Donaldson theory is 
gapless [Witt5j . Moreover, both the N = 2 and N = 4 supersymemtric (3 + 1)- 
TQFTs are not unitary. It seems to be an open question whether or not there are 
unitary and gapped 3 + 1-TQFTs that can detect smooth structures of 4-manifolds. 

In an exposition for mathematicians |Witt6j . Witten speculated about a possible 
connection between Seiberg- Witten theory and superconductors. The challenge 
to better understand (3 + l)-TQFTs lies at the frontier in both topology and the 
exploration of topological phases of matter. A realization of Seiberg- Witten theory 
in condensed matter systems would be a landmark in the interaction of topology 
and physics. 
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